Controlled excitation and resonant acceleration of ultracold few-boson systems by 

driven interactions in a harmonic trap 
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We investigate the excitation properties of finite utracold bosonic systems in a one-dimensional 
harmonic trap with a time-dependent interaction strength. The driving of the interatomic coupling 
induces excitations of the relative motion exclusively with specific and controllable contributions of 
momentarily excited many-body states. Mechanisms for selective excitation to few-body analogues 
of collective modes and acceleration occur in the vicinity of resonances. We study via the few- 
body spectrum and a Floquet analysis the excitation mechanisms, and the corresponding impact 
of the driving frequency and strength as well as the initial correlation of the bosonic state. The 
fundamental case of two atoms is analyzed in detail and forms a key ingredient for the bottom- 
up understanding of cases with higher atom numbers, thereby examining finite-size corrections to 
macroscopic collective modes of oscillation. 

PACS numbers: 67.85.-d,05.30. Jp,03.75.Kk,03.65.Ge 



I. INTRODUCTION 

Detailed control of both single particle potential land- 
scapes and interparticle interactions is an appealing fea- 
ture of ultracold atom physics ^ . Apart from cooling the 
atoms to ultra low temperatures, it is in the meanwhile 
routinely possible to design almost arbitrarily shaped op- 
tical and magnetic traps and to tune them with an un- 
precedented control [2:]. The interactions among ultra- 
cold atoms can be adjusted by exploiting magnetically 
and optically induced Feshbach resonances Q. The di- 
mensionality of the system can be tuned by strongly con- 
fining the transverse degrees of freedom leading to quasi 
one-dimensional traps and dimensionality-specific phe- 
nomena like the Tonks Gas [1, Q of impenetrable bosons 
and its attractive excited state counterpart, the Super 
Tonks gas J,]. For these low dimensional systems an ad- 
ditional tool to control the interactions are confinement 
induced resonances @. 

In the present work we take advantage of this exper- 
imental progress concerning the design of external and 
most importantly interatomic forces, and consider the ef- 
fects of a time-dependent oscillating interaction strength 
in a one-dimensional harmonic trap. Time-dependent 
driving is usually applied to external traps, providing in- 
spiring effects such as the dynamical control of tunneling 
[13 1 dynamical localization flTl , photon-assisted tun- 
neling [13] via a periodic driving of the lattices or the 
excitation of collective oscillations in a harmonic trap 
[isj to mention only a few. However, investigations con- 
sidering a time-dependent scattering length, usually re- 
ferred to as 'Feshbach resonance management' [Tj] for 
the mean-field situation, have inspired a lot of research 
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on control of solitons [l5| or modulational instabilities 
[l6l |. Experimental investigations in this direction have 
been performed recently [13, [ll]. The main advantage 
of the driving of the scattering length compared to other 
driving modes applied on the external potential for exam- 
ining collective excitations [l3, il^ > that other species 
or non-condensate fractions are not affected by the driv- 
ing of the interaction of one-species [3] ■ Apart from the 
harmonic trap, a two mode system with time-varying in- 
teraction has been studied with Floquet theory, leading 
to many-body coherent destruction of tunneling and lo- 
calization (20| . 

On the other hand, beyond the mean-field regime, 
there are relatively few works dealing with the time- 
dependent modulation of the scattering length, address- 
ing mainly the experimental results on the formation 
of molecules [13] from a few-body perspective [2l|, |2^ . 
While these works concentrate mainly on the attractive 
part of the two-body spectrum and on the correspond- 
ing bound state, our work focuses exclusively on repul- 
sive interactions. As illustrated in Ref. |2l] in a har- 
monic trap the coupling of an excited two-body state to 
the molecular ground state is very efficient since nearby 
states are out of resonance, and possess a relatively small 
coupling to the initial state. For the repulsive case as we 
will demonstrate, the relevant few-body states reflect the 
equidistant spectrum of the harmonic trap and therefore 
lead to a more complex dynamics involving several in- 
stantaneous configurations. A recent publication (23j ex- 
plored the integrability of the system, via a similar model 
with time-modulated interaction thereby calculating the 
dynamical structure factor. 

In this work, we aim to examine from a few body 
perspective the effects resulting from a periodic modu- 
lation of the repulsive interaction strength, focusing on 
few-body collective excitations, control of the dynam- 
ics and state population, as well as mechanisms of ac- 
celeration via resonances. The case of two atoms in 
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an one-dimensional harmonic trap for which the energy 
spectrum is known analytically for a constant arbitrary 
strength of the interaction [23] serves as a starting point 
of the investigation. The successful experimental prepa- 
ration of few-body systems with a controllable number 
of atoms has been realized in Ref. [l^ and the ener- 
gies have been measured with high precision (26j . Since 
the modulation of the interaction strength affects exclu- 
sively the relative motion in this system (the center of 
mass is decoupled and therefore unaffected), we study 
explicitly the internal motion. Firstly the focus is on the 
frequencies and the driving amplitudes that give rise to 
a controlable excitation to particular states after prepa- 
ration in a certain initial state with a specific strength 
of the interaction. Additionally we examine particular 
few-body analogues of collective macroscopic modes of 
breathing oscillations, as well as a resonant acceleration 
mechanism via multiple excitations. Our analysis of res- 
onances and acceleration modes is supported by calcula- 
tions of the Floquet spectrum for the effective single de- 
gree of freedom of the relative motion within a harmonic 
trap and an oscillating delta barrier. Going to higher 
atom numbers we demonstrate similarities and analyze 
the differences with the basic case of two particles, and 
compare the results with those of macroscopic calcula- 
tions, showing finite size effects on the collective modes. 
All calculations are performed by the numerically exact 
Multi-Configurational Time-Dependent Hartree method 
(MCTDH see Appendix) , which is especially designed to 
treat the dynamics of many degrees of freedom under 
time-dependent modulations. 

This article is organized as follows: In Section II we in- 
troduce our model, and in Section III our focus is on the 
case of two particles thereby examining the mechanisms 
of controlable collective excitations to specific states and 
the infiuence of parameter changes on them. We investi- 
gate in Section IV the acceleration mechanism via mul- 
tiple excitations and calculate the Floquet spectrum of 
this case illustrating the underlying mechanism for the 
appearance of the resonances. An extension of this study 
to higher atom numbers is performed in Section V con- 
centrating on the analogue of the breathing mode for 
collective oscillations and finite size corrections. In the 
last Section VI we summarize our results and provide an 
outlook. 



II. MODELING THE TIME DEPENDENT 
INTERACTION 

Quasi one-dimensional waveguides can be created 
by choosing a strongly focused laser field yielding 
strongly confined transversal directions compared to 
the longitudinal one. In this way the trap becomes 
highly anisotropic with the characteristic length for the 

transversal trapping a± e 



much smaller than the 



longitudinal one (a;j_ is the frequency of the harmonic 
confinement in the transversal direction). Consequently 



the transverse degrees of freedom are energetically frozen 
as only the ground state is occupied and the effective ID 
interaction strength reads for the case of contact interac- 
tions 3: 
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where the free-space s-wave scattering length oq does not 
depend on the detailed appearance of the potential for 
the interatomic interaction, which is then modeled by an 
effective contact potential. There are two parameters in 
EqUthat can be tuned to attain a time-dependent inter- 
action strength: (i) the scattering length ao via a change 
of the strength of e.g. a magnetic field B approaching to 
or departing from a Feshbach resonance -Feshbach res- 
onance management- as oq = o,bgi^ — ) '^ti^re Ai? 
and Bq are the width and the position of t6e resonance, 
respectively, and ahg is the background scattering, and 
(ii) the transversal length a± by modifying the relevant 
laser parameters, taking into account the quasi-one di- 
mensional restrictions posed above (see also Ref. [isl ). 

The one-dimensional A^-body Hamiltonian with a time 
dependent coupling g(t) reads: 
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where we have performed a scaling transformation set- 
ting the length scale equal to the longitudinal charac- 
teristic oscillator length a\\ = ^/jj^ and the energy 



scale to /iwjj, while the scaled interaction strength is 
g — 9id\/mJWuj\^. The interaction potential between 
each pair of particles i,j is represented by the Dirac S- 
function. We note that for numerical purposes, a Gaus- 
sian with a very small width of the order of the grid 
spacing is employed. 

Initially, the particles are prepared in the ground state 
of the harmonic trap with an interaction strength . We 
will then explore the excitation dynamics for a periodic 
driving of the repulsive interaction strength of the form: 



ait) = .90 + Agsin^iujt) 



(3) 



where Ag is the amplitude of the driving and uj the driv- 
ing frequency. The impact of each of these three pa- 
rameters of the driving law (gQ,Ag,uj) will be examined. 
The reason for the specific choice sin^ (ujt) for the driving 
is our focus on repulsive interactions, i.e., g should stay 
positive even for 170 = 0. Since sin^(a;) — L-£2^^2l±L^ the 
above driving law comprises a periodic oscillation with 
frequency 2u}. Inve stig ating purely attractive interac- 
tions as done in Ref. [2l| or alternating between attractive 
and repulsive interactions such as in Refs. [3, [H rep- 
resent interesting but different situations. 
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III. RELATIVE MOTION OF THE TWO ATOM 
PROBLEM AND INSTANTANEOUS 
EIGENSPECTRUM 

In general, N particles in a harmonic trap with contact 
interaction represents a separable problem H = Hj^ + Hr- 

N 

The center of mass R — jj^ Xi and and its conjugate 

i=l 

^ d 

momentum P = —i} — — constitute the center of mass 

i—l 

Hamiltonian iJ/j = ^ + . The Hamiltonian of the 
relative motion Hr is in general not subject to further 
simplifications and cannot be solved analytically. Never- 
theless, for the special case of two particles the relative 
motion (r = xi — X2) reduces to an effective one-body 
problem: 

Hr=pl + 'j+gmr) (4) 

The contact interaction affects only the relative motion, 
leaving unaffected the center of mass. Therefore we fo- 
cus on the relative motion, which actually represents a 
one particle problem with a harmonic trap and a delta 
barrier with oscillating height placed in the center. With 
the transformation r — ^/2x, g{t) = V2g'{t) we obtain 
the standar^ form of the harmonic oscillator Hamilto- 
nian H'j, = ^ + ^+g'{t)S{x). Hr defines an analytically 
solvable eigenvalue equation [23| in the case of a time- 
independent parameter g, which we will discuss next as 
it is very important for the understanding of the excita- 
tion dynamics of the relative motion. The solutions cover 
in general the complete interval g € [0, 00) and we will 
refer to them as the instantaneous eigenstates where 
n = 0, 1, 2... are the energy levels for a certain time in- 
stant tref- In spite of the existence of these stationary 
solutions the driven time-dependent problem possesses 
no closed analytical solution, although a study via the 
evolution of the coefficients in an expansion with respect 
to the corresponding instantaneous eigenstates is natural 
^21, 22]. 

In Fig. [ija) we show the lowest lying eigenenergies for 
Hr ■ For g — the eigenspectrum of Hr is the harmonic 
oscillator single-particle spectrum with the equidistant 
eigenenergies e„ = 71 + 1/2. As g increases the odd levels 
are unaffected by S{x) since they possess a node at the co- 
ordinate origin, while the even states acquire an increas- 
ing energy and a dip at x = 0. Therefore, as we observe 
in Fig. [TJa) each even level approaches energetically the 
next upper odd level forming a doublet spectrum, char- 
acteristic for double-well potentials. In the limit where 
the 'barrier height' is infinite g ^ 00 the even levels be- 
come degenerate with the odd ones. This limit g — >■ 00 
is the so-called Tonks- Girardeau limit where the bosons 
are mapped to non- interacting fermions ■ The Tonks- 
Girardeau Gas, is one of the most fundamental systems 
appearing exclusively in one-dimensional many-body sys- 
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Figure 1: (a) The three lowest eigenenergies of the symmetric 
(0,2,4) and antisymmetric (1,3,5) eigenstates of the relative 
motion for two particles with increasing interaction strength 
g. (b) The energy difference between the eigenstates of the 
relative motion £2 — eo and £4 — £2 with increasing g. 

tems. The transition to 'fermionization' for two atoms 
has been recently reported also experimentally [26l |. en- 
hancing the interest on few-body studies, like the present 
one. 

The energy difference between two even levels with in- 
creasing interaction strength shown in Fig. [IJb) plays 
also a crucial role for the dynamics. In fact, since the ini- 
tial preparation is in the ground state which corresponds 
to an even state of the relative motion, the dynamics 
can only lead to a population of other even states, and 
the corresponding energy distance is crucial for the time- 
evolution. For the two limits of zero interaction and the 
Tonks-Girardeau gas (g — > 00) the gap between two even 
parity levels is |e„ — e„_2|c(=o , g-j-oo = 2. Starting from 
g — and increasing g this value slightly decreases [see 
Fig. [ijb)], since the states with larger quantum number 
n possess a lower probability density at a; = and are 
therefore less affected by the contact interaction. The 
slope at g = reads: 

^LH0.Wx)I«H*.(O)P^;;^^ (5) 

The response to a minor increase of g therefore depends 
on the value of the harmonic oscillator eigenstates at 
a; = which decreases with n [see Fig. H^a)]. Therefore 
at the onset of the interactions the successive even states 
tend to approach each other energetically. Fig. [Ijb) 
shows the effect of an increasing coupling strength g on 
the distance between the first two even levels eo and €2 
(red line). The most rapid change of this energy gap 
is near g — with a minimum value £2 — eo ~ 1-85 at 
g* w 2.2 and it approaches asymptotically the value 2 for 
(7 — >■ 00. The slope ^^^^ is decreasing as g increases [see 
Fig. lUa)] and from some value of g on, the energy gap 
starts to increase again and asymptotically approaches 
the value of the non-interacting system. Additionally, 
since decreases with increasing n the deviation from 
the value 2 for £2 — eo shown in Fig. [Ijb) (red line) is 
the largest possible such deviation between two succes- 
sive (even) states. This is exemplarily shown in the same 
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figure by the distance between the next pair of succes- 
sive states 62 and 64 (green fine). It is obvious that the 
energy gap 64 — £2 is always larger than the £2 — eo gap, 
and this holds analogously also for gaps between higher 
lying neighboring states. 

The above discussed features of the energy spectrum 
and the corresponding gaps, possess a crucial impact on 
the dynamics which we will discuss later. As g{t) changes 
with time according to Eq. ^ different regions of g- 
values in Fig. [T] are probed according to the choice of 
the parameters go and Ag. We might already foresee e.g. 
that a driving around small g-values possesses a greater 
impact than for a driving for larger 17- values, since the 
corresponding slope is larger, subsequently leading to 
larger energy variations. The equidistance of the spec- 
trum close to the two extreme limits, as well the decrease 
of the gaps at small to intermediate values of g will also 
be of great importance concerning a possible resonant 
behaviour (see sections IV and V). 

We note that Fig.[T]is based on the exact eigenenergies 
obtained from the equation [23 |: 
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The corresponding even eigenstates are given analytically 
in the form of parabolic cylinder functions. For the nu- 
merical calculations of the time-dependent evolution of 
the system we use a regularized delta-function of the 



form: Sa{x) = 



with a = 0.05 which is small 



enough to catch the real 'delta-like' behavior but also 
convenient for a numerical grid sampling. Minor numeri- 
cal deviations on the eigenstates and the quantum evolu- 
tion stemming from this approximation are unavoidable. 



IV. RESONANT CONTROLLABLE 
EXCITATION 

The wave packet of the relative motion at t — cor- 
responds to the ground state of Hr in Eq. (jH) with 
g{0) — ga- This an even state and since parity is con- 
served during the time evolution only even parity states 
can be occupied, which relates to the bosonic permuta- 
tion symmetry. In this work we are interested mainly in 
two quantities: 

• the population of momentarily eigenstates 4>n at a 
certain time p„ (ire/) = |(0n(ire/)|'/'(^re/))P, where 
tp{t) is the wave function of the relative motion. We 
use tref = (fc € No) without loss of generality. 

• the time evolution of the energy e{t). We refer only 
to the relative energy since the center of mass is 
completely untouched by the change of g. 

We start our investigation in the regime of small driv- 
ing amplitudes where the excitation dynamics is to a 



Figure 2: Population of instantaneous eigenstates with time 
variation, for (a) go = 0.0, Ag = 0.05, uj = 0.99, (b) go = 0.0, 
Ag = 0.05, iJ = 0.997, (c) go = 2.0, Ag = 0.1, cj = 0.925. 
Profile of the resonances via the minimal occupation of the 
ground state with varying driving frequency. Several cases 
are shown for different values of go and Ag. 



larger extent controllable. The main role in this regime is 
played by the driving frequency. If this frequency is much 
lower than the gap between two even states (w << 1) 
then we are in the adiabatic regime and the evolution 
involves only the momentarily ground state (j)Q (with 
Po > 0.99). Approaching the first resonance lo = 1 from 
below induces an excitation to the state (j>2 of the rela- 
tive motion. This is shown for a typical case in Fig. [2] (a) 
where go = 0, Ag = 0.05 and oj = 0.99. We see that the 
ground state looses population while the second excited 
state gains. The next even level 04 remains almost un- 
populated. For the particular case though where go = 
(and correspondingly for go — > oo) since the eigenspec- 
trum is initially completely equidistant, an excitation to 
the next level (/)4 is not prohibited corresponding to a two 
step process. Therefore we see in Fig [3] (b) that a bit 
closer to the resonance (w = 0.997) the n = 4 level gains 
population after the n — 2 level does so. This is also in 
general the case for larger amplitudes Ag > 0.5 where 
multiple excitations are enhanced as we will see later on. 
Therefore, while the system departs ~ in this case com- 
pletely - from the ground state, the first collective excita- 
tion to state n = 2 is necessarily combined with a transfer 
of population to the next level and therefore a control- 
lable excitation exclusively to the second level (complete 
depopulation of ground state and complete population of 
the second excited level) is not possible here. It becomes 
though possible if the initial interaction is stronger. For 
instance for go = 2.0, Ag = 0.1 we see in Fig. ^ (c) a 
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complete transfer of population to the first excited level. 
Therefore, the non-harmonicity in the spectrum due to 
the initial correlations is helpful from the point of view 
of a controllable state excitation and preparation. Let us 
note that for a controlable creation of such states, one 
should choose a small amplitude < 0.5, since larger 
amplitudes lead easily to multiple occupation of excited 
states due to the close to equidistant spectrum. Addi- 
tionally the driving frequency should be carefully tuned 
to be close to the corresponding energy spacing of the 
spectrum with 17(0) = Qq. For this case uj — 0.925 - 
somewhat lower than the resonant frequency for the non- 
interacting go = situation- we have a decreased gap in 
the energy spectrum [see Fig. [T] (b)]. 

As an overview of the resonances in this regime we 
present in Fig. [2] (d) the minimal occupancy of the in- 
stantaneous ground state po as a function of the driving 
frequency for several small amplitudes of the driving. We 
observe that far from resonance the minimal po ~ Ij so 
there is hardly any excitation as expected. The frequency 
where the ground state becomes at a certain time com- 
pletely unoccupied corresponds to the resonance. The 
resonant frequency is slightly shifted to lower values for 
a larger amplitude of the driving which is attributed to 
the decrease of the energy gap in the spectrum [see Fig. [T] 
(b)], as Ag covers larger regions of g (for go = 0). We 
also observe that for go = 2.0 the resonant frequency is 
shifted to much lower values than for go = 0, approxi- 
mately corresponding to the energy gap of the levels at 
g = 2.0 (w = w 0.925). For this case though, small 

changes in the driving amplitude do not shift significantly 
the position of the resonance, since the energy spacing 
in the vicinity of this value of go does not change sub- 
stantially [see Fig. [T](b)]. However, the most important 
difference of the two cases for different initial coupling 
strength go is the fact that the case go = leads easily 
to multiple excitation of higher states since the energy 
gaps are all close to the same value corresponding to the 
same resonant frequency, while for larger go the energy 
spacing is not equidistant and therefore a complete con- 
trolable transfer to a certain state is possible. 

A comment on the robustness of the initial state prepa- 
ration is in order here. Apart from being easily ex- 
citable to different instantaneous states, the initially non- 
interacting ensemble is in general more sensitive to the 
driving of the interaction. Even far from resonance the 
evolution of this initial state {go = 0), leads to a change 
in energy of the order of 10% while for stronger and par- 
ticular intermediate interactions as well as close to the 
fermionization regime go = 20.0 it is ten orders of mag- 
nitude lower. This is understandable if one inspects the 
slopes of the energy curves in Fig [T] (a), which are much 
larger for small values of g. This could be a signature for 
the detection of a highly correlated ensembles like the 
Tonks Girardeau gas, i.e., by studying their response to 
changes of the interaction strength. 

Not only an excitation to the first excited state of 
the relative motion is possible, but also to other excited 
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Figure 3: Occupation of the adiabatic eigenstates of the 
Hamiltonian for go — 0, Ag — 0.5 close to the second and 
the first resonance (a) uj = 1.98 (b) u — 0.998 

states if the resonant frequency is chosen correspondingly 
as we can see in Fig. |3] (a). As expected the resonance 
width though decreases for transitions to higher states. 
The effect of the initial go and of the amplitude Ag is 
similar to the previous case. 

For larger amplitudes the controllability of the excita- 
tion process reduces, as many states are subsequently ex- 
cited, and simultaneously taking part in a complex time 
evolution. A typical example is presented in Fig. [3] (b). 
Still the frequency plays the dominant role and only close 
to resonances the evolution leads to highly excited states 
of the spectrum. Mechanisms of acceleration appear then 
which we will discuss in the following section. We would 
like to note here that a large Ag, offer the possibility 
of 'multi-photon excitations'. In this case even for low 
frequencies which are an integer ratio of the principle 
resonant frequency, excitations become possible. 



V. ACCELERATION VIA MULTIPLE 
EXCITATION AND FLOQUET ANALYSIS 

We will now more thoroughly examine the case of 
strong driving, which makes it possible, as we will show, 
to accelerate the particles, i.e., to increase the mean value 
of the energy with time. The process of multiple excita- 
tions as we have seen, is possible close to resonances, since 
the spectrum is approximately equidistant. Especially a 
larger value of Ag leads to a covering of wider areas of the 
energy gaps, and therefore the comparatively small differ- 
ences between the gaps effectively drain away. Through 
this multiple excitation process, the system never returns 
completely to the ground state, and indeed occupies grad- 
ually increasingly higher lying states. This excitation 
process induces an increase of the energy to very high 
values as long as the gaps to higher excited states are in 
the resonance window. 

We present in Fig. 2] the time evolution of the expec- 
tation value of the energy, close to the first resonance 
for Ag — 20.0. For values of the frequency sufficiently 
far from resonance, a repopulation of the ground state 
in the course of the time can be observed, while a multi- 



6 



50 



40 



0) 
\- 

0) 

c 
LU 



30 



20 



(0=1.05 
(0=1.0 

(0=0.95 
(0=0.9 
(0=0.8 



20 40 60 80 100 

time 




Figure 4: Time evolution of the expectation value of the en- 
ergy for Ag = 20.0, go — and several frequencies close 
to resonance. Acceleration shown close at the resonant fre- 
quency 



mode behavior is encountered due to an excitation of 
several states. Approaching the resonance, the instan- 
taneous ground state is never repopulated significantly, 
on the opposite, higher states of the spectrum are subse- 
quently populated in the same manner as shown in Fig. 
[3] (b). This leads to an acceleration, i.e., energy gain of 
the particles. Our finite time simulations indicate that 
this energy gain approaches a saturation to very high 
values of the energy for strong but finite Ag. 

Let us analyze this resonant mechanism from the per- 
spective of Floquet theory which has been developed for 
time-periodic Hamiltonians. The Ansatz of the Floquet 
theory for the time-dependent wave packet reads: 



'>Pk{r,t) 



(7) 



where u''{r, t) = J^n it^(?')e*"'^* are the so-called Floquet 
eigenstates or quasi-energy states, which are time peri- 
odic functions expanded here into a Fourier series and qk 
are the quasienergies. Introducing this Ansatz into the 
Schrodinger equation for the Hamiltonian in Eq. Q we 
find: 

gfc<(r) = Inuj + + \r^) -f (go + \Ag)6{T)\ ui{r) 



(8) 



The r.h.s. of this equation is the Floquet Hamiltonian 
which has diagonal, i.e., proportional to ujij, and off- 
diagonal terms coupling the different modes u^_2, 'U'n+2- 
To find ujj and eigenvalues qk one should solve this eigen- 
value problem by diagonalizing the corresponding Hamil- 
tonian. The harmonic oscillator fm{i~) basis is very con- 
venient to express the Floquet Hamiltonian in our par- 
ticular problem, not only because of the harmonic part 
of the potential which has as matrix elements the har- 
monic oscillator eigenvalues but also because the matrix 




Figure 5: Floquet eigenspectrum for (a) increasing driving 
frequency u with go — 0.5 Ag = 0.5, (b) increasing initial 
interaction strength go for the resonant frequency lo — 1 and 
Ag = 0.5 



elements of 5{r) are of the form J drf,n{r)S{r)f„i' (r) = 
/m(0)/m'(0). An alternative method equally well-suited 
for our problem is to write the Floquet eigenvectors in 
terms of parabolic cylinder functions [27| D{—a — \tX) 
which are the solutions of the underlying stationary prob- 
lem (23 | . In this representation we have 



with the important properties D{a, 0) 



(9) 
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By integrating Eq. JHl us- 

ing these relations we obtain a recursive system of alge- 
braic equations for the coefficients c„. Demanding that 
this system possesses solutions we obtain the values of 
quasienergies qk- We have used both methods to derive 
the eigenspectrum of the Floquet Hamiltonian for our 
problem and confirmed the agreement crosschecking the 
results and thereby the numerical convergence. 

We show in Fig. [S] (a) the eigenspectrum of the Flo- 
quet Hamiltonian, i.e., the Floquet quasi-energies with 
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Figure 6: The lowest eigenenergies with increasing interaction 
g for (a) three and (b) four bosons. 

increasing driving frequency lo for a typical case go = 0-0 
and A (7 — 1.0. In the dense spectrum of quasi-energies we 
encounter points of accumulation when the frequency is 
at resonant values (w = 1,2,3...). At these points, many 
Floquet quasi-energies take a value close to the harmonic 
oscillator eigenenergies (0.5, 1.5, 2.5...), and form Floquet 
bands. The even- frequency resonances (w = 2,4,6...) 
correspond to accumulation points with an energy gap 2, 
i.e. at the quasi-energy values 0.5,2.5,4.5, .... 

Let us now show how these results obtained for the 
Floquet Hamiltonian are connected with the quantum 
acceleration processes for quantum resonances. In dy- 
namical systems a pure point Floquet spectrum is as- 
sociated with localized behavior and the energy remains 
bounded at all times, while singularly continuous compo- 
nents are responsible for diffusive behavior and growth of 
the energy (see Ref. ^28*] and references therein). At the 
resonant driving frequencies as we have seen above the 
quasi-energies are accumulating and approach particular 
values forming close to continuous areas. This property 
of the spectrum leads to the acceleration and energy gain. 

Besides, in Fig. [5] (b) we see that with increasing 50 
and for the resonant frequency ui = 1 the eigenenergies 
of the Floquet spectrum deviate from each other, and 
only come close again as we approach the fermioniza- 
tion limit in the next upper level. This deviation from 
the accumulation points makes the evolution less diffu- 
sive for intermediate interactions, and corresponds to the 
picture of the instantaneous spectrum, where the energy 
gaps for intermediate values of g become less equidistant, 
prohibiting multiple excitation. 



VI. HIGHER ATOM NUMBERS AND FINITE 
SIZE EFFECTS ON COLLECTIVE EXCITATIONS 
(BREATHING MODE) 

The extensive analysis of the two-body case above, is 
very useful for an understanding of the effects occurring 
for higher atom numbers. The main reason for this are 
the properties of the many body spectrum of the har- 
monic oscillator including the delta-type interaction. We 
present in Fig. [5] the energetically low lying part of the 
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Figure 7: (a) population of instantaneous eigenstates for 
three particles with go = 2.0, Ag = 0.05 and u) — 0.9. 
(b) Population of the instantaneous ground state for differ- 
ent numbers of particles N = 2,3,4 go — 0.0 Ag = LO and 
Lu = 0.6 



spectrum for three and four particles with increasing in- 
teraction strength [see also Ref. [lil . We observe that all 
states show a quite similar evolution and thus the energy 
gaps between them do not deviate significantly (cross- 
ings or anti-crossings do not occur). Also important is 
that most of the states correspond to excitations of the 
center of mass, which are not relevant to our study. For 
example the first excited state and one of the two states 
of the second excited band which behave exactly like the 
ground state, correspond to the ground state of the rela- 
tive motion and to the first and second excitation of the 
center of mass motion, respectively. 

The time-dependent variation of the interaction 
strength which affects only the relative motion offers the 
possibility of a controllable excitation to specific states 
also for higher atom numbers, since states like those men- 
tioned above which correspond to excitations of the cen- 
ter of mass do not contribute to the time evolution and 
are therefore avoided. We demonstrate this in Fig. [7] 
(a) for the case of three particles and an initial value 
50 = 2.0 in the intermediate regime, where an excitation 
predominantly to the lower state (corresponding to an 
excitation of the relative motion) of the second excited 
band of eigenstates is performed [see Fig. EJa)]. We de- 
note here by the numbers 0, 2, 4 the energetically ordered 
states which correspond to a respective excitation of the 
many body relative motion. The resonant frequency for 
excitation from state to state 2 will be referred as the 
principle resonance in the following. 

The above signifies that higher number of particles al- 
low for a similar control of their dynamics as in the case 
for two particles which is mainly due to the similarities 
of their underlying energy spectrum. One important dif- 
ference here is that for equal parameters, systems with 
larger atom numbers experience a stronger impact of the 
driving than those with a lower atom number. This can 
be seen from the the maximum loss of population of the 
instantaneous ground state with time in Fig. [7] (b) in 
a non-resonant case which gets larger with an increas- 
ing number of particles. An explanation for this atom 
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number related effect is the response of the ground state 
energy to the variation of g: 

dE N{N ~ 1) /", , 

where we see that the slope of the total energy at g = 
increases quadratically with the particle number. There- 
fore a variation of g possesses a greater impact for higher 
atom numbers. 

Another important observation concerning the size of 
the system is the position of the principle resonance in 
the regime of intermediate interactions. We have seen 
already that for e.g. go = 2.0 the resonant frequency is 
lower than for very weakly or strongly interacting initial 
states. This frequency becomes even lower as the num- 
ber of particles increase which is based on the decreas- 
ing energetical spacing in the corresponding many-body 
spectra with increasing particle number. Consequently a 
lower value of frequency is needed for larger atom num- 
ber for the corresponding resonant excitation via driving 
of the interaction strength. 

The above observation is important for modes of col- 
lective oscillation of the wave function, in analogy to the 
macroscopic collective oscillations f?]. Many measure- 
ments in experiments are based on exciting collective 
modes, which are usuall y a nalyzed within the effective 
mean field descriptions [I^. The frequencies of these 
oscillations can be obtained with a high accuracy from 
the experimental data usually observing the size or the 
mean position of the condensate 0, [13 j and represent a 
very important measure for identifying different interac- 
tion regimes. For example in one dimension the Thomas 
Fermi limit, or the Tonks Girardeau and Super- Tonks gi- 
rardeau gas, possess a very characteristic ratio between 
the so-called dipole mode which is an oscillation of the 
center of mass of the condensate with the trap frequency 
idd, and the breathing or first compression mode, where 
essentially the size of the condensate is oscillating with 
a frequency Wb. For the effectively non-interacting limits 
{g = 0,oo) this ratio is ujb/ujd — 2 while for the Thomas 
Fermi hmit it is Wb/w^ = ^3 « 1.73 p^. These im- 
portant theoretical results have been confirmed in the 
corresponding experiments 0, [13 ■ 

From the few body perspective which we examine here, 
the mode of oscillation that we excite by varying the 
scattering length, is of a compressional-breathing type. 
We demonstrate this in Fig. |8] where we show several 
snapshots of the one-body density for characteristic time 
moments. The excitation to the second state of the spec- 
trum which corresponds to an excitation of the relative 
motion (and therefore to a broader wave function), pos- 
sesses therefore the characteristics of a breathing mode. 
The characteristic frequency of such a mode, in the in- 
termediate interaction regime could be compared with 
the mean field result for the Thomas Fermi regime. For 
non- interacting cases we confirm the ratio Wfc/wt = 2. 
The Thomas Fermi regime applies to large ensembles 
of particles, and the analogy with a the finite system 
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Figure 8: One body density for the case of three particles as 
in Fig. [7] for several snapshots. 

examined here can be at most indicative. Neverthe- 
less in one-dimensional systems we can take the min- 
imum value of the energy gap which appears close to 
g = 2.0 as corresponding to the Thomas-Fermi regime 
(see also corresponding arguments in [soj). In doing 
so, we can reveal finite size corrections for the breath- 
ing mode of the Thomas Fermi limit, starting with two 
particles Wb/wt ~ 1.85 [Slj and going to Wb/wt ~ 1.78 for 
five particles with a further tendency to decrease with 
increasing particle number. In the macroscopic limit 
— ^ cx) we shoul d ap proach the mean field estimate 
WbM = V3 w 1.73 [l9l]. 

This interrelation of the collective mode frequencies 
and their finite size deviations within the few-body spec- 
trum are probably valuable for further studies on con- 
necting, checking and reinterpreting mean field results 
in the light of the exact many-body spectrum and be- 
haviour. 



VII. CONCLUSION 

We examined the effects of a periodically driven in- 
teraction strength on ultracold bosonic systems in a one- 
dimensional harmonic trap which can be realized by time 
a modulation of magnetic fields utilizing Feshbach reso- 
nances or periodically changing the transversal confine- 
ment length in a waveguide via laser fields. We have 
shown that the feature of near equidistant energy levels 
for arbitrary interaction strength for the relative motion 
of two particles in a harmonic trap but also for the corre- 
sponding many-body spectrum for larger atom numbers 
has important consequences for the excitation dynam- 
ics. In particular, the energetical spacing yields resonant 
driving frequencies, which one can employ to excite par- 
ticular states of bosonic relative motion. We underline 
that unlike other driven many-body systems, the varia- 
tion of the interaction strength offers the possibility to de- 
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sign excitations of tlie relative motion of a certain species 
exclusively with a very high degree of controllability for 
certain regimes of the driving parameters. Approaching 
the resonant frequencies the atoms excite to the corre- 
sponding excited level of the relative motion which has 
been demonstrated here by calculating the occupation of 
instantaneous eigenstates. For strong driving amplitudes 
in the vicinity of the resonances the energy reaches out 
to very high values with several and successive excita- 
tions to energetically higher lying states of the spectrum. 
This multi-excitation process of acceleration is also an- 
alyzed via the properties of the Floquet spectrum. The 
initial interaction strength distorts the energy spectrum 
shifting the position of the resonances, while highly cor- 
related initial states are quite insensitive with respect to 
the changes of the repulsive interaction. We have shown 
via our exact numerical calculations, that for any num- 
ber of particles there is a similar and for larger ensembles 
even more sensitive response to the driving of the inter- 
action strength leading to higher excitation amplitudes. 
Effects due to the finite size of the system, are also ana- 
lyzed from the perspective of collective oscillation modes, 
and especially the analog to the macroscopic breathing 
mode is established, thereby discussing similarities and 
deviations from mean field approaches. Especially the 
two-body problem discussed here represents a case of in- 
terest to the experiment of distinguishable fermions in a 
harmonic trap |25l [26j . Interesting outlooks are the ex- 
ploration of different driving laws, including the possibil- 
ity to alter between repulsive and attractive interactions 
or using different potential landscapes. 
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Appendix A: Computational Method 

Treating time-dependent Hamiltonians with many in- 
teracting degrees of freedom is a computationally very 
demanding problem. In this work for all numerical cal- 
culations we rely on the Multi-Configurational Time- 
Dependent Hartree (MCTDH) method [32. - .34J , primarily 
a wave-packet dynamical tool known for its outstand- 
ing efficiency in high-dimensional applications. The un- 
derlying idea of MCTDH is to solve the time-dependent 
Schrodinger equation 



as an initial-value problem by an expansion in terms of 
direct (or Hartree)products $j: 



using a convenient multi-index notation for the configu- 
rations, J — {ji . . . jf), where f — N denotes the number 
of degrees of freedom and Q = (xi , . . . , Xf)'^ . The single- 
particle functions ip^J^^ are in turn represented in a fixed, 
primitive basis implemented on a grid. For indistinguish- 
able particles as in our case, the single-particle functions 
for each degree of freedom k = 1 , . . . , are of course 
identical in both type and number {(pj^, with < n). 

In the above expansion, both the coefficients Aj and 
the Hartree products are time-dependent. Using the 
Dirac-Frenkel variational principle, one can derive equa- 
tions of motion for both Aj, $ j. This conceptual compli- 
cation offers an enormous advantage: the basis t)} 
is variationally optimal at each time t, allowing us to keep 
it fairly small. Exactly for this reason MCTDH is ide- 
ally designed to solve time-dependent Hamiltonians like 
driven systems as the one we tackle here. The permu- 
tation symmetry can be enforced by symmetrizing the 
coefficients Aj, but the ground state is automatically of 
bosonic character. 

In addition the Heidelberg MCTDH package [35] incor- 
porates the so-called relaxation method which provides a 
way to obtain the lowest eigenstates of the system by 
propagating some wave function ^Pq by the non-unitary 
g--f/r (^pj-QpdgdiiQji imaginary time.) As r — cxd, this 
automatically damps out any contribution but that stem- 
ming from the true ground state. 



-H 



-EjT 



|J)(J|*o) 



i^ = 
^'(Q,0)-vi/o(g) 



(Al) 



In practice, one relies on a more sophisticated scheme 
termed improved relaxation [36]. Here {^\H — E\^i) is 
minimised with respect to both the coefficients Aj and 
the configurations $j. The equations of motion are 
solved iteratively, first for Aj{t) (by diagonalisation of 
($j]i7]<i>if) with fixed $j) and then propagating $j in 
imaginary time over a short period. The cycle will then 
be repeated. This method is used in the present work 
to compute the eigenstates and spectrum for A = 3, 4, 5 
particles. 

As it stands, the effort of this method scales exponen- 
tially with the number of degrees of freedom, . Just 
as an illustration, using 15 orbitals and N = 5 requires 
7.6 • 10^ configurations J. This restricts our analysis in 
the current setup to about N = 0(10), depending on 
how decisive correlation effects are. Therefore we are 
limited in terms of numerical convergence to the num- 
ber of particles we can treat for large variations of the 
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interaction strength. The number of orbitals needed for 
a time-evolution of our driven system also depends on 
how prominent the excitations are, i.e., how close to res- 
onance the driving frequency is. The latter fact sets lim- 
itations to examining the long-time dynamics where the 
behaviour involves many excitations rather than oscilla- 
tory. By contrast, the dependence on the primitive ba- 



sis, and thus on the grid points, is not as severe. In our 
case, the grid spacing should of course be small enough 
to sample the interaction potential. The length of the 
grid should be also choosen sufficiently large especially in 
cases where the driving frequency approaches resonances 
and the time-dependent wave function covers a very wide 
coordinate range. 
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